Shell model analysis of competing contributions to the double-beta decay of 48 Ca 
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Background: Neutrinoless double beta decay, il observed, would reveal physics beyond the Stan- 
dard Model (SM) of particle physics, namely it would prove that neutrinos are Majorana fermions 
and that the lepton number is not conserved. 

Purpose: The analysis of the results of neutrinoless double beta decay observations requires an 
accurate knowledge of several nuclear matrix elements (NME) for different mechanism that may 
contribute to the decay. We provide a complete analysis of these NME for the decay of the ground 
state (g.s.) of 48 Ca to the g.s. 0^ and first excited 0^ state of 48 Ti. 

Method: For the analysis we used the nuclear shell model with effective two-body interactions that 
were fine-tuned to describe the low-energy spectroscopy of p/-shell nuclei. We checked our model 
by calculating the two-neutrino transition probability to the g.s. of 48 Ti. We also make predictions 
for the transition to the first excited 0% state of 48 Ti. 

Results: We present results for all NME relevant for the neutrinoless transitions to the 0^ and 0^ 
states, and using the lower experimental limit for the g.s. to g.s. half-life we extract upper limits 
for the neutrino physics parameters. 

Conclusions: We provide accurate NME for the two-neutrino and neutrinoless double beta decay 
transitions in A=48 system, which can be further used to analyze the experimental results of double 
beta decay experiments when they become available. 

PACS numbers: 23.40.Bw, 21.60.Cs, 23.40.-s, 14.60.Pq 

Keywords: Double beta decay, Nuclear matrix elements, Shell model 
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I. INTRODUCTION 

Neutrinoless double beta (0vj3(3) decay, which can only 
occur by violating the conservation of the total lepton 
number, if observed it will reveal physics beyond the 
Standard Model, and it will represent a major milestone 
in the study of the fundamental properties of neutri- 
nos [l|-[7|. Indeed, its discovery would decide if neu- 
trinos are their own antiparticles [8|, and would provide 
a hint about the scale of their absolute masses. That 
is why there are intensive investigations of this process, 
both theoretical and experimental. Recent results from 
neutrino oscillation experiments have demonstrated that 
neutrinos have mass and they can mix However, 
the neutrino oscillations experiments cannot be used to 
determine the neutrino mass hierarchy and the lowest 
neutrino mass. Neutrinoless double beta decay is viewed 
as one of the best routes to decide these unknowns. A key 
ingredient for extracting the absolute neutrino masses 
from Ov/3/3 decay experiments is a precise knowledge of 
the nuclear matrix elements (NME) for this process. 

There are potentially many mechanisms that could 
contribute to the neutrinoless double beta decay pro- 
cess that will be briefly reviewed below. Several of these 
mechanisms do not provide contributions to the decay 
rate that explicitly depend on the neutrino masses, but 
their effect would vanish if the neutrinos are not massive 
Majorana particles Q. In all cases the half-lives depend 
on the nuclear matrix elements that need to be accurately 



calculated using low-energy nuclear structure models. In 
particular, if the exchange of light left-handed neutri- 
nos is proven to be the dominant mechanism, one could 
be able to use the experimental results and the associ- 
ated NME to extract the neutrino mass hierarchy and 
the lowest neutrino mass Q. The two- neutrino double 
beta (2v/3f3) decay is an associate process that is allowed 
by the Standard Model, and it was observed in about 
ten isotopes. Therefore, a good but not sufficient test of 
nuclear structure models would be a reliable description 
of the 2v(3fi half-lives. 

Since most of the j3/3 decay emitters are open shell nu- 
clei, many calculations of the NME have been performed 
within the pnQRPA approach and its extensions (T^|-[23j. 
However, the pnQRPA calculations of the more common 
two-neutrino double beta decay half-lives, which were 
measured for about 10 cases 24), are very sensitive to the 
variation of the so called g pp parameter (the strength of 
the particle-particle interactions in the 1 + channel) [121 ]- 
[l4| , and this drawback still persists in spite of various 
improvements brought by its extensions EH -[111, includ- 



ing higher-order QRPA approaches [2l|-[23|. The out 
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come of these attempts was that the calculations became 
more stable against g pp variation, but at present there 
are still large differences between the values of the NME 
calculated with different QRPA-based methods, which 
do not yet provide a reliable determination of the two- 
neutrino double beta decay half-life. Therefore, although 
the QRPA methods do not seem to be suited to predict 
the 2i/j3j3 decay half-lives, one can use the measured 2u/3f3 
decay half-lives to calibrate the g pp parameters, which 
are further used to calculate the Of/3/? decay NME [2a |. 
Other methods that were recently used to provide NME 
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for most Ov/3/3 decay cases of interest are the Interacting 
Boson Model (IBM-2)j2p,IIz!], the Projected Hatree-Fock 
Bogoliubov (PHFB) [2J and the Generator Coordinate 
Method (GCM) [29|. 

Recent progress in computer power, numerical algo- 
rithms, and improved nucleon-nucleon effective interac- 
tions, made possible large scale shell model calculations 
(LSSM) of the 2v(3/3 and Qi/(3j3 decay NME [30|-[32|. The 
main advantage of the large scale shell model calcula- 
tions is that they seem to be less dependent on the ef- 
fective interaction used, as far as these interactions are 
consistent with the general spectroscopy of the nuclei in- 
volved in the decay. Their main drawback is the limi- 
tation imposed by the exploding shell model dimensions 
on the size of the valence spaces that can be used. The 
most important success of the large scale shell model cal- 
culations was the correct prediction of the 2vj3/3 decay 
half-life for 48 Ca [3(J HH]. In addition, these calcula- 
tions did not have to adjust any additional parameter, 
i.e. given the effective interaction and the Gamow- Teller 
(GT) quenching factor extracted from the overall spec- 
troscopy in the mass-region (including beta decay proba- 
bilities and charge-exchange strength functions) , one can 
reliably predict the 21/0/3 decay half- life of 48 Ca. 

Clearly, there is a need to further check and refine these 
calculations, and to provide more details on the analysis 
of the NME that could be validated by experiments. We 
have recently revisited [H| the 2v(3/3 decay of 48 Ca using 
two recently proposed effective interactions for this mass 
region, GXPF1 and GXPF1A, calculating the NME and 
half-lives for the transition of the 48 Ca g.s. to the g.s. 
and the first excited 2 + state of 48 Ti. 

In this paper we add to the analysis the 2i/f3/3 transi- 
tion to the first excited state of 48 Ti. We also extend 
our analysis (3|| of the Ov/3/3 decay of 48 Ca by provid- 
ing the NME associated with the most important 0i/(3/3 
mechanisms for transitions to the g.s. and first excited 
O2 state of 48 Ti. Future experiments on double beta de- 
cay of 48 Ca (CANDLES j37| and CARVEL [H}) may 
reach the required sensitivity of measuring such transi- 
tions, and our results could be also useful for planning 
these experiments. 



II. TWO-NEUTRINO DOUBLE BETA DECAY 

LSSM calculations of 2i/{3(3 decay NME can now be 
carried out rather accurately for many nuclei [Hj] . In the 
case of 48 Ca, Ref. [3(| reported for the first time a full pf- 
shell calculation of the NME for the 2vf3/3 decay mode, 
for both transitions to the g.s. and to the 2^ excited 
state of 48 Ti. As an effective interaction it was used the 
Kuo-Brown G-matrix [40j with minimal monopole mod- 
ifications, KB3 41 1 . In Ref. [13] we use the recently 
proposed GXPF1A two-body effective interaction, which 
has been successfully tested for the pf shell 42|-[4J|, to 
perform 2v(3f3 decay calculations for 48 Ca. Our goal was 
to obtain the values of NME for this decay mode, for both 



transitions to the g.s. and to the 2+ state of 48 Ti, with 
increased degree of confidence, which would allow us to 
consider similar calculations for the Ou/3/3 decay mode of 
this nucleus [Hj]. The 2i//3/3 transitions to excited states 
have longer half-lives, as compared with the transitions 
to the g.s., due to the reduced values of the correspond- 
ing phase space factors, but they were measured in some 
cases, such as 100 Mo |45| . 

For the 2v/3f3 decay mode the relevant NME are of 
Gamow- Teller type, and has the following expression for 
decays to states in the grand-daughter that have the an- 
gular momentum J = 0, 2 [l|-@, 



1 



(^ + lkr-||l+)<l+||<7r-||0+) 



VJ + T 



(E k +Ej)J+i 



(1) 

Here E^ is the excitation energy of the it state of inter- 
mediate odd-odd nucleus, and Ej — ^Qpp(J + ) + AM. 
Qp@{J + ) is the Q- value corresponding to the 0/3 decay 
to the final J~j state of the grand-daughter nucleus, and 
AM is the mass difference between the parent and the 
intermediate nucleus 48 Sc. The most common case is the 
decay to the 0^ g.s. of the grand-daughter, but decays to 
the first excited 0^ and 2\ states are also investigated. 
The 2i/(3/3 decay half-life expression is given by 



T, 



2u,J 

1/2 



G¥\M*? T (J)\ 2 



(2) 



where G 2 j are 2vj3ft phase space factors. Specific values 
of G 2 f for different 2v/3f3 decay cases can be found in 
different reviews, such as Ref. 



of G 2 r see Ref. [46|. In Ref. [3 



For a recent analysis 
j occ j-uoj.. puij . ±11 ira. [o-aj we explicitly analyzed 
the dependence of the double-Gamow- Teller sum enter- 
ing the NME Eq. ([1]) vs the excitation energy of the 1 + 
states in the intermediate nucleus 48 Sc. This sum was re- 
cently investigated experimentally (35j , and it was shown 
that indeed, the incoherent sum (using only absolute val- 
ues of the Gamow- Teller matrix elements) would provide 
an incorrect NME, thus validating our prediction. We 
have also corrected by several orders of magnitude the 
probability of transition of the g.s. of 48 Ca to the first 
excited 2+ state of 48 Ti reported in Ref. [3(j. 

In Ref. 34j we fully diagonalized 250 1 + states in the 
intermediate nucleus to calculate the 2v/3ft decay NME 
for 48 Ca. This procedure can be used for somewhat heav- 
ier nuclei using the J-scheme shell model code NuShellX 
[HI, but for cases with large dimension one needs an al- 
ternative method. The pioneering work on 48 Ca [3C| used 
a strength-function approach that converges after a small 
number of Lanczos iterations, but it requires large scale 
shell model diagonalizations when one wants to check the 
convergence. Ref. 49] proposed an alternative method, 
which converges very quickly, but it did not provide a 
complete recipes for all its ingredients, and it was never 
used in practical calculations. Recently [13], we proposed 
a simple numerical scheme to calculate all coefficients of 
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TABLE I: Matrix elements and half-lives for 2u decay calcu- 
lated using GXPF1A interaction and two quenching factors. 
Matrix elements are in MeV -1 for transitions to + states 
and in MeV -3 for transitions to 2 + states. 





qf 


= 0.77 


qf 


= 0.74 


TTT 


M 2v 


i% (y) 


M 2v 


T?h (y) 


ot 


0.054 


3.3 x 10 19 


0.050 


3.9 x 10 19 




0.012 


8.5 x 10 23 


0.010 


1.0 x 10 24 


Ot 


0.050 


1.6 x 10 24 


0.043 


1.9 x 10 24 



the expansion proposed in Ref. [4!| . Following Ref. [4!| , 
we choose as a starting Lanczos vector, Lf, either the 
initial or final state in the decay (only + to + tran- 
sitions are considered), to which we apply the Gamow- 
Teller operator. This approach is very efficient for large 
model spaces, as for example the j j55 space (consisting of 
the 0g7/2, Id, 2s, and h\\i-x orbits), which for the 128 Te 
decay leads to m-scheme dimensions of the order of 10 
billions necessary to calculate the g.s. of 128 Xe. In the 
calculation of 48 Ca decay we use the standard quenching 
factor, qf — 0.77, for the Gamow- Teller operator ot. We 
checked the result reported in Ref. [34| using this alter- 
native method and we found the same result. The novel 
result report here for the first time is for the transition 
to the first excited 0+ state in 48 Ti at 2.997 MeV. The 
matrix element when using GXPF1A interaction is 0.050, 
very close to that for the transition to the g.s. Using the 
phase space factor Gf + = 2.43 x 10~ 22 MeV 1 from Ref. 

Q (a new set of phase space factors were recently pro- 
posed [46|, but for 2vj3fi decays they differ only by 4% 
from those of Ref. Q) we found that the half-life for this 
transition is 1.6 x 10 24 y. We recall here that our results 
reported in [34| for the half-lives of the transitions to g.s. 
and to the first 2 + excited state are 3.3 x 10 19 y and 
8.5 x 10 23 y, respectively. One can see that the transition 
to the first excited Oj" state at 2.997 MeV is predicted to 
compete with the transition to the first excited 2^~ state 
at 0.994 MeV. 

The half-life for the transition to the g.s. 0^~ was mea- 
sured by several groups with increased precision (see e.g. 
[24j]). The most recent result from NEMO-3 (see |24| and 
references therein) is Ty 2 = AA^o'listat.) ± 0A(syst.). 
Our GXPF1A result is marginally out of the recently 
reduced error bars. However, a recent publication [50] 
found a quenching factor of 0.74 for the p/-shell nuclei 
using GXPF1A interaction. The same quenching factor 
was proposed some time ago [5lJ using a different ef- 
fective interaction. Using the smaller quenching factor 
of 0.74 brings the calculated half-life within the experi- 
mental limits. A comparison of the matrix elements and 
the associated half-lives for the two quenching factors 



used here is given in Table I. Potential observation of the 
21/(9/3 transitions to the excited states of 48 Ti could shed 
some light on the variation of the quenching factor for 
the Gamow- Teller operator in this nucleus. One should 
also mention that the excitation energy of the 0^~ state 
in 48 Ti calculated with GXPF1A interaction is about 1 
MeV higher than the experimental value, while it is about 
right for 48 Ca. Other available effective interactions do 
no provide a better description of this state. This re- 
sult may raise concerns about the validity of the nuclear 
structure description of this state within the p/-shell. An 
experimental observation of the 2vf3(3 transition to this 
state could be used to validate (or not) our result. 



III. NEUTRINOLESS DOUBLE BETA DECAY 

The Ov/3/3 decay, (Z, A) -> (Z + 2, A) + 2e~, requires 
the neutrino to be a massive Majorana fermion, i.e. it 
is identical to the antineutrino [8j. We already know 
from the neutrino oscillation experiments that some of 
the neutrinos participating in the weak interaction have 
mass, and that the mass eigenstates are mixed by the 
PNMS matrix Uik, where I is the lepton flavor and k is 
the mass eigenstate number (see e.g. Ref. |52j). How- 
ever, the neutrino oscillations experiments cannot decide 
the mass hierarchy, the mass of the lightest neutrino, and 
some of the CP non-conserving phases of the PNMS ma- 
trix (assuming that neutrinos are Majorana particles). 

Considering only contributions from the exchange of 
light, left-handed(chirality), Majorana neutrinos Q, the 
Ov(3f3 decay half-live is given by 

T^r 1 = G 0iy |M°f (i^MX) 2 . (3 ) 

J \ Tfl e J 

Here, G 0l/ is the phase space factor, which depends on 
the 0v(3f3 decay energy, Qpp, the charge of the decaying 
nucleus Z, and the nuclear radius [3, S(|. The effective 
neutrino mass, (mpp), is related to the neutrino mass 
eigenstates, rrik, via the left-handed lepton mixing ma- 
trix, U ek , 

{mpp) /m e = ri vL = ^ m kUl k /m e . (4) 

k— light 

m e is the electron mass. The NME, M® v , is given by 

M„°" = M(jt - (j^) M° F U - M° v , (5) 

where M^ T , Mf and Mfi are the Gamow-Teller (GT), 
Fermi (F) and tensor (T) matrix elements, respectively. 
Using closure approximation these matrix elements are 
defined as follows: 
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M°J = (0+ \Y,r- m r- n 0^ n | 0+ \ 

\ m,n / 

J] TBTD (j pJp , , JnJn , ; r) (j pJp , ■ rT I r_ 1 r_ 2 0« I j„j n , ; J^T) a , (6) 

jp3p>3nj n lJ-K 



where 0" m are Oz//3/3 transition operators, a = 
(GT, F, T), | 0+ > is the g.s. of the parent nucleus, 
and | 0y > is the final + state of the grand daugh- 
ter nucleus. The two-body transition densities (TBTD) 
can be obtained from LSSM calculations 36]. Expres- 
sions for the anti-symmetrized two-body matrix elements 
(TBME) (j P jpr,J*T\ T - 1 T- 2 0%\'U n r,J*T) a can be 
found elsewhere, e.g. Refs. [3g, [53[. Assuming that 
one can unambiguously measures a Qu(3j3 half-life, and 
one can reliably calculate the NME for that nucleus, one 
could use Eqs. and (UJ) to extract information about 
the lightest neutrino mass and the neutrino mass hierar- 
chy [521 ] . In addition, one could consider the contribution 
from the right-handed currents to the effective Hamilto- 
nian, which can mix light and heavy neutrinos of both 
chiralities (L/R) 



where Nk are the heavy neutrinos that are predicted by 
several see-saw mechanisms for neutrino masses [52| . Uik 
and Vik are the left and right-handed components of the 
unitary matrix that diagonalizes the neutrino mass ma- 
trix [54j . One should also mention that there are several 
other mechanisms that could contribute to the 0v/3f3 de- 
cay, such as the exchange of supersymmetric (SUSY) par- 
ticles (e.g. gluino and squark exchange [55j]), etc, whose 
effects are not directly related to the neutrino masses, 
but indirectly via the Schechter-Valle theorem Q. As- 
suming that the masses of the light neutrinos are smaller 
than 1 MeV and the masses of the heavy neutrinos, Mk, 
are larger than 1 GeV, the particle physics and nuclear 
structure parts get separated, and the inverse half-life 
can be written as 



VeL = ^ U ^kL + ^ U ekN k L 
k— light k— heavy 

VeR = ^ V ^ v kR + ^ V ^ N kR , (7) 
k— light k— heavy 



1 l/2 



= G°" \r) uL M^ + < A > X\+ < rj > X n + ( VNL + VNR ) M°" 



VkkM° k » k \ 2 . 



(8) 



where t] v l was defined in Eq. (|4|), and 

EjtI m p 
ek M~k' 

k— heavy 

( M Wl \ A ^ 2 m p 
VNR * \M^J „ f 

k— heavy 

< A > = £ ^ U ekVek, 
k— light 

x ' k— light 

Here e is the mixing parameter for the right heavy bo- 
son Wr and the standard left-handed heavy boson Wl, 
Wr « eWi + W2, Mwr and Mwl are their respective 



masses, and m p is the proton mass. The rjy and r/g are 
the R-parity violation contributions in supersymmetric 
(SUSY) Grand Unified Theories (GUT) related to the 
long range gluino exchange and squark-neutrino mech- 
anism, respectively [52] . Finally, the t}kk term is due 
to possible Kaluza-Klein (KK) neutrino exchange in an 
extra-dimensional model [5a]. The set of nuclear ma- 
trix elements M°", X\, X v , , My, and M\ v are 
discussed in many reviews, e.g. Ref. (52[. The M^ K 
analysis can be found in Ref. (56|. In particular, using 
the factorization ansatz [56| one gets 

r,KKM° K \ = <m>SA AC + ™ P < ™ X > M° N » 

= VikkAC + r lhKK M° N », (10) 
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where < m >sa and < m~ x > KK masses depend on the 
brane shift and bulk radius parameters, and are given in 
Table II of [56]. One can see that the mass parameters 
< m >sa /fn e and m p < mT 1 > has the effect of modify- 
ing and r)NR respectively. | m p < mT x >|< 1CU 8 and 
it could in principle compete with t]mr- |< m >sa /m e | 
varies significantly with the model parameters and it 
could also compete with t) v l- One needs to go beyond 
the factorization ansatz, and use information from several 
nuclei (57j to discern any significant contribution from the 
KK mechanism. 

Constraints from colliders experiments suggest that 
terms proportional with the mixing angles, e, f7 e fe(7ieauy) > 
and V e k(ught) are very small [54j. The present limits are 
|< A >|< 1CT 8 and |< rj >|< 1CT 9 , but they are expected 



to be smaller. In addition, the contributions from X\ and 
X v terms in Eq. ([8]) would produce angular and energy 
distribution of the outgoing electrons different than that 
coming from all other terms (H , and these signals are un- 
der investigation at SuperNEMO [581 ] . Here we assume 
that these contributions are small and can be neglected. 
In addition, if < A > is small, Eq. ((9|) suggests that 
t)nl is small. Information from colliders also puts some 
limits on (M Wr , M n ) ~ (2.5GeV, lAGeV), and these 
limits will be refined at LHC in the coming years. Based 
on this information and the present limit on the Oi//3/3 
decay of 76 Ge one can estimate that | rj V L |< 10~ 6 , and 
I Vnr |< 10~ 8 . Then, the half-life can be written as 



J l/2 



rOv 



rQvl 



(11) 



r 



TABLE II: Matrix elements for Ov decay using GXPF1A in- 
teraction and two SRC models [6lj], CD-Bonn (SRC1) and 
Argonne (SRC2). For comparison, the (a) values are taken 
from Ref. [27(. and the (b) value is taken from Ref. [Kg] for 
g pp = 1 and no SRC. 







M° N V 


M° X Y 


Mf 




SRCl 


0.90 


75.5 


618 


86.7 




SRC2 


0.82 


52.9 


453 


81.8 




others 


2 g(a) 


46.3 (a) 


392 (b) 






SRCl 


0.80 


57.2 


486 


84.2 




SRC2 


0.75 


40.6 


357 


80.6 



where we adjusted r\ v h and r\jqR for potential KK contri- 
butions, f] vL = r/„ L + rii KK and ij N = r/ NR + r\ hKK . 

If one neglects the SUSY and KK contributions until 
a hint of their existence is provided by colliders experi- 
ments or future results of 0v(3/3 decay experiments show 
that these contributions are necessary [57| . then 



x; 



1/2 



= G " ( \M? V 



2 | ,2 



M 



A'* 



\Vnr\ 



(12) 

where we used the fact that the interference between the 
left-handed terms and the right-handed terms is negligi- 
ble [HI. 

The structure of the M^ 1 is the same as that described 



H a (r) (see e.g. page 68 of Ref. [Hj]). A detailed descrip- 
tion of the matrix elements of Of 2 for the jj-coupling 
scheme consistent with the conventions used by modern 
shell model effective interactions is given in Ref. [36| . 
One should also mention that our method [36[ of calcu- 
lating the TBTD. Eq. (|6j), is different from that used in 
other shell model calculations [32|]. We included in the 
calculations the recently proposed higher order terms of 
the nucleon currents, three old and recent parametriza- 
tion of the short-range correlations (SRC) effects, finite 
size (FS) effects, intermediate states energy effects, and 
we treated careful few other parameters entering the into 
the calculations. We found very small variation of the 
NME with the average energy of the intermediate states, 
and FS cutoff parameters, and moderate variation vs the 
effective interaction and SRC parametrization. We could 
also show that if the ground state wave functions of the 
initial and final nucleus can be accurately described us- 
ing only the valence space orbitals, the contribution from 
the core orbitals can be neglected. This situation is dif- 
ferent from that of the nuclear parity-nonconservation 
matrix elements (59|, for which the "mean- field" type 
contribution from the core orbitals could be significant 
[60j |. Another important result that clearly transpires 
from our formalism is that in the closure approximation 
the neutrinoless transition to the first excited 2 + state 
is zero. This result is due to the rotational invariance 
of the TBME entering Eq. © (see also Appendix of 
Ref. [36j]). The structure of the R- parity breaking SUSY 
mechanisms NME is similar to that of light and heavy 
neutrino exchange mechanisms, but with no a = F com- 
ponent [Hj|. The neutrino potentials used here for the 
Mjy", and those used for the most significant contribu- 
tions to M» and M~ v NME are given in Ref. [H, but 



in Eqs. (5)-(8), with slightly different neutrino potentials for completeness they are reviewed in the Appendix with 
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TABLE III: Single mechanism upper limits for neutrino 
physics parameters rjj extracted from the lower limit of the 
half-life for the transition to the ground state of 48 Ti [S^] and 
using the matrix elements from Table II. 





\Vul\ x 10 5 


\W\ x 10 7 


M X 10 8 


\ m \ x 10 7 




SRC1 


3.79 


4.52 


5.52 


3.93 




SRC2 


4.16 


6.45 


7.53 


4.17 



the specific parameters included in these calculations. 

The results for all NME entering Eq. (JTTJ) for the tran- 
sition to the g.s. and first excited 0^" state of 48 Ti are 
presented in Table II. Comparison with results of other 
models, when available, are also included. For the light 
neutrino exchange matrix element we choose to compare 
with the IBA-2 results, which is very different from ours. 
Other shell model analyses of this particular NME gives 
similar results for both transitions to Of and 0^ states 
[iU, |63|. To our knowledge, with the exception of the light 
neutrino exchange NME, no other results of shell model 
calculations for these matrix elements were reported so 
far (with the possible exception of Ref. [64{ where the 
NME as a function of neutrino mass is reported and it 
could potentially be used to extract the corresponding 
Mjtf). Based on these calculations and using the exper- 
imental lower limit of the half-life, one can extract the 
" single- mechanism dominance" upper limits for | rjj |, 
where j = {vL), N, A', q. At present there is available 
only the lower limit of the half-life for the transition to the 
g.s. of 48 Ti 1 4x 10 22 y [52j|. Using the phase-space factor 
from Ref. [46], G 0u = 61.4 x 1CT 15 y 1 (for g A = 1.254 
and R = 1.2A 1 / 3 fm), we obtained the upper limits for 
| rjj | shown in Table III. Alternatively, assuming that 
two or more mechanisms are contributing to the half-life 
in Eq. (|lip compete, one could use the experimental data 
from several isotopes to assess the contribution of each 
mechanism [55| . Clearly, this scenario requires as many 
as possible accurate half-lives and associated NMEs. For 
example, in the likely scenario that no more than two 
mechanisms are competing and they are the light and 
heavy neutrino exchange, then Eq. (fl"2"|) can be used to 
analyze the data. If the exchange of light neutrino will be 
determined as the dominant mechanism, then our results 
could possible be used to decide the light neutrino mass 
hierarchy and the lowest neutrino mass |52j |. 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we analyzed the 2^/3/3 and several mech- 
anisms that could compete to the 0^/3/3 decays of 48 Ca 
using shell model techniques. We described very efficient 
techniques to calculate accurate 2vj3/3 NME for cases that 



involve large shell model dimensions. These techniques 
were tested for the case of 48 Ca, and we provided NME 
and half-lives for 2z//3/3 transitions to the g.s. and ex- 
cited states of 48 Ti. These techniques can be used to 
make predictions for 76 Ge, 82 Se using the jj44 model 
space (O/g/2, lp, 0g 9/2 ), and for 128 Te, 130 Te and 136 Xe 
using the model space. 

We reviewed the main contributing mechanisms to the 
OvPP decay, and we showed that based on the present 
constraints from colliders one could reduce the contribu- 
tion to the 0i/j3j3 half-life to the relevant terms described 
in Eq. (fTTj) . A reliable analysis of the Qvf5f3 decay exper- 
imental data requires accurate calculations of the associ- 
ated NME. We extended our recent analysis [3(| of the 
0^/3/3 NME for 48 Ca to include the heavy neutrino ex- 
change NME, the long range gluino exchange NME, and 
the squark-neutrino mechanism NME. We also presented 
for the first time shell model results of these new NME 
for the 0^/3/3 transitions to the g.s. and the first excited 
0^ state in 48 Ti. 

To extend this analysis to the A > 48 cases, more ef- 
forts have to be done to include all spin-orbit partners in 
the valence space and satisfy the Ikeda sum-rule, reduce 
the center-of-mass spurious contributions, and better un- 
derstand the changes in the effective 0^/3/3 transition op- 
erators [HI, [66| • In addition, the closure approximation 
used to calculate the NME within the shell model ap- 
proach and by other methods (e.g. IBA-2 [26[, PHFB 
[HI, and GCM |29j) needs to be further checked for accu- 
racy, especially for the heavy neutrino exchange, the long 
range gluino exchange, and the squark-neutrino mecha- 
nism. An analysis of the double beta decay of 136 Xe that 
addresses some of these issues is in preparation. 



V. APPENDIX 

The matrix elements for the light and heavy neutrino 
exchange in Eq. (fTTj) have the same structure as that 
described in Eqs. (3)-(6) of Ref. [H- For M ?/ the 
neutrino potential is the same as in Eq. (7) of [36j 

H a (r) = — ^ f a (qr)MQ-G a ( q *) q dq, (13) 
ir J q+ (E) 

with the same ingredients described in Eqs. (9)-(12) of 
[36j . Here we corrected the (fx p — fj, n ) value to 4.71, an 
error that seems to be propagating for some time through 
the literature Q. This correction explains the small dif- 
ference between the M® v values of Table II and corre- 
sponding ones reported in Ref. [36j. Fortunately, this 
correction only changes the matrix elements by few per - 
cents. All other constants are the same as in Ref. |36| . 
In particular, we used g A = 1.254 and R = 1.2 A 1 / 3 fm. 
For the there is a slight change in the neutrino po- 
tentials 
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H a (r) 



2R 



irm e m p J 



fa(qr)h a (q 2 )G a (q 2 )q 2 dq,(U) 



where m e and m p are the electron and proton mass, re- 
spectively. 

The most significant contributions to M®Y and M~ u 
have a similar structure as and M^, however, only 
the a = GT, T terms in Eq. ([5]) are contributing. The 
radial neutrino potentials for M^Y have the same form as 
those used for Mjy", Eq. ([13]), but with different h a : 



The radial neutrino potentials for M~ v have the same 
form as those used for M°", Eq. (IT3"1) . but with different 
h„: 



1 

T-GT,T — 



q 2 /ml 



6 m e (m„ + m d ) (1 + q 2 jm\Y 



(16) 



where m n and are the current up and down quark 
masses. In the calculation we used m u + = 11.6 
MeV. 



l GT,T 



m e m p q 2 /m* 2m e m p q 2 /m* 



l + q 2 /ml (l + q yml) 2 

(15) 

where m T is the charged pion mass, 139 MeV. Expres- 
sions for c l7r and c 2jr are given in Ref. The numerical 1068217 and the SciDAC Grant NUCLEI is acknowl- 
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